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A  Proposed  Gas  Release  Experiment 
on  the  ARGOS  Satellite 


1.  INTRODUCTION 


Artificial  modification  and  ionization  enhancement  in  the  ionosphere  may  result  in  local 
ionospheric  irregularities  that  may  control  or  interfere  with  radio  wave  propagation,  generate 
tell-tale  infrared  signatures,  enhance  or  degrade  existing  ionospheric  irregularities,  and 
produce  a  wide  variety  of  phenomena  through  non-linear  processes.1  In  our  proposed 
experiment,  we  plan  to  release  xenon  and  carbon  dioxide  gases  into  the  atmosphere  from  the 
Advanced  Research  and  Global  Observation  Satellite  (ARGOS)  to  study  (1)  the  critical 
ionization  velocity  (CIV)  process2  of  xenon  gas  and  (2)  the  chemical  reactions  of  0O2  gas.  In 
this  report,  we  discuss  the  main  points  of  the  proposed  experiment  and  present  a  feasibility 
study. 


Received  for  publication  24  January  1992 

1  Bernhardt.  P..  Tepley.  C.A.,  and  Duncan,  L.M.  (1989)  Airglow  enhancements  associated  with 
plasma  cavities  formed  during  ionospheric  heating  experiments,  J.  Geophys.  Res.,  94:9071. 

2  Alfven,  H.  (1960)  Collision  between  a  non-ionlzed  gas  and  a  magnetized  plasma.  Rev.  Mod. 
Phys.,  32:710-713. 
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2.  THE  ARGOS  SATELLITE 


The  satellite  will  be  in  a  circular  orbit  at  an  altitude  of  about  450  nmi.  (835  km)  with  an 
orbital  velocity  of  7.4  km/s.  The  gas  release  experiment  will  be  one  of  several  experiments  to 
be  conducted  onboard  the  satellite.  It  will  be  conducted  when  the  satellite  is  passing  over 
Maui,  Hawaii,  where  the  AMOS  telescope  will  observe  the  experiment  (Figure  1). 


3.  XENON  GAS  RELEASE 


Xenon  gas  will  be  released  in  four  different  directions  (Figure  2).  The  purpose  of  the 
experiment  is  to  study  the  signatures  of  enhanced  ionization  due  to  the  CIV  process. 


3.1  Critical  Ionization  Velocity 

Alfven3  proposed  the  hypothesis  of  critical  ionization  velocity  (CIV)  -  when  a  neutral  gas 
travels  through  a  magnetized  plasma  with  a  relative  velocity  exceeding  a  critical  value,  rapid 
ionization  of  the  neutral  gas  occurs.  The  critical  value  of  the  velocity  is  called  the  critical 
ionization  velocity  V*.  and  is  given  by:3 


V.  =  V(2e4>/M) 


where  e4>  is  the  ionization  energy  and  M  is  the  mass  of  a  neutral  particle  (Figure  3).  V.  is  the 
minimum  value  of  Vj_  to  cause  ionization,  where  Vx  is  the  relative  velocity  component 
perpendicular  to  the  magnetic  field.  For  xenon,  e«t>  =  12.13  eV,  M  =  131.3  AMU.  and  V.  = 
4.225  km/s. 


3  Alfven,  H.  (1954)  On  the  Origin  of  the  Solar  System,  Oxford  Univ.  Press,  Oxford. 

/ 


NETIC 


Figure  1 .  The  Proposed  ARGOS  Satellite  CIV  Experiment  Over  Maui 
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Figure  2.  Directions  of  Gas  Releases  From  ARGOS  in  the  Proposed  Experiment 
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Figure  3.  A  Schematic  Diagram  of  CIV 
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Moblus  et  al4  first  suggested  that  since  spacecraft  velocity  may  be  close  to  V.,  the  vector 
sum  of  the  spacecraft  and  gas  velocities  may  equal  or  exceed  V,.  Axnas5  worked  out  some 
necessary  conditions  for  a  CIV  Interaction  between  the  ionospheric  plasma  and  a  xenon  cloud. 
Murad  et  al6  suggested  releasing  four  gases  from  the  Space  Shuttle  for  testing  CIV  In  space. 


3.2  Plasma  Instability 

The  CIV  process  Is  generally  accepted  to  be  mediated  by  some  beam-plasma  Instability, 
such  as  the  two-stream  plasma  Instability7  8  9  that  occurs  when  a  neutral  beam  (with  Ions 
and  electrons  In  It)  travels  across  an  ambient  magnetic  field  (Figure  4).  The  beam-plasma 
Instability  waves,  energized  by  the  ion  beam,  accelerate  electrons.  According  to  quasi-linear 
plasma  theory,  the  electrons  form  a  plateau-tall  distribution.10  Some  of  the  'hot'  electrons  In 
the  tall  may  be  energetic  enough  to  Ionize  the  neutrals  by  impact.  As  Ionization  occurs,  new 
Ions  are  created.  When  they  are  created,  they  have  the  velocity  of  the  neutrals.  They  carry  on 
the  CIV  process  by  energizing  the  beam-plasma  Instability  waves,  which  In  turn  energize  the 
electrons.  Thus,  CIV  Is  a  cyclic  process,  in  which  a  neutral  gas  is  Ionized  by  means  of  Its  own 
kinetic  energy. 


4  Moblus,  E.,  Boswell,  R.W.,  Piel,  A.,  and  Henry,  P.  (1979)  A  spacelab  experiment  on  the 
critical  ionization  velocity,  Geophys.  Res.  Lett.,  6:29. 

5  Axnas,  I.  (1980)  Some  necessary  conditions  for  a  critical  velocity  interaction  between  the 
Ionospheric  plasma  and  a  xenon  cloud,  Geophys.  Res.  Lett,  7:933-936. 

6  Murad,  E.,  Lai,  S.T.,  and  Stair,  A.T.  (1986)  An  experiment  to  study  the  critical  ionization 
velocity  theory  in  space,  J.  Geophys.  Res.,  91:A9  10188-10192. 

7  Moblus,  E.,  Papadopoulos,  K.,  and  Plel,  A.  (1987)  On  the  turbulent  heating  and  the 
threshold  condition  In  the  critical  ionization  velocity  Interaction,  Planet.  Space  Science, 
35:345. 

8  Goertz,  C.K.,  Machida,  S.,  and  Lu,  G.  (1990)  On  the  theory  of  CIV,  Adv.  Space  Res.,  33-46. 

9  McNeil.  W.J.,  Lai,  S.T.,  and  Murad,  E.  (1990)  Interplay  between  collective  and  collisional 
processes  in  critical  velocity  Ionization,  J.  Geophys.  Res.,  95:A7  103435-10354. 

10  Formlsano,  V.,  Galeev,  A.A.,  and  Sagdeev,  R.Z.  (1982)  The  role  of  the  critical  ionization 
velocity  phenomenon  in  the  production  of  the  Inner  coma  cometary  plasma.  Planet.  Space  Set, 
30:491-497. 
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Table  1 .  Important  Natural  Parameters  for  the  CIV  Experiment  on  ARGOS 


Xenon  Mass 

131.3  AMU 

Xenon  Ionization  Potential 

12.13  eV 

Xenon  Critical  Ionization  Velocity  V* 

4.225  km/s 

Altitude 

835  km 

Nighttime  Ambient  H+  Density 

2X104  /cc. 

Daytime  Ambient  0+  density 

2X104  /cc. 

Nighttime  Ambient  Electron  Density 

2X104  /cc. 

Daytime  Ambient  Electron  Density 

5x10®  /cc. 

Magnetic  Field  Magnitude  at  Maui 

0.36  G 

Magnetic  Field  Declination 

11° 

Magnetic  Field  Inclination 

38.7° 

Electron  Gyrofrequency 

1006x10®  Hz 

Xenon  Gyrofrequency 

4.17  Hz 

Xenon  Gyroradius 

282.4  m 

Xenon  Lower  Hybrid  Frequency 

2.047  kHz 

3.3  Seed  Ionization 

To  begin  with,  the  neutral  gas  must  have  seed  ionization,  otherwise  no  plasma  interaction 
will  occur.  In  the  proposed  experiment,  there  is  no  thermal  ionization  by  explosions,  stripping 
ionization  is  energetically  impossible  at  the  gas  velocity  used,  and  photoionization  is  ruled  out 
because  the  ionization  potential  (12.13  eV)  of  xenon  is  higher  than  the  Lyman-a  line  energy 
(10.2  eV)  of  sunlight.  Seed  ions  can  come  only  from  charge  exchange  with  the  ambient  plasma. 
In  the  daytime,  the  charge  exchange1 1  is  between  Xe  and  ambient  0+,  the  cross  section  being  o0+ 
=  4.8xl0~16  cm2.  The  density  of  0+  at  835  km  altitude  is  typically  2xl04/cc.  At  night,  the 
charge  exchange  is  between  Xe  and  ambient  H+,  the  cross  section12  being  oH+  =  1.5xl0-16  cm2 
and  density  of  H+  =  2xl04/cc.  Assuming  that  the  ions  leave  the  beam  in  1/4  gyro  period  tb,  the 
number  density  of  xenon  ions  generated  by  charge  exchange  is  given  by 

rV4 

ni(t)=  dtnxe°Vnamblent  (2) 

Jo 


11  Jonsson,  B.O.  (1980)  private  communication  to  I.  Axnas. 

12  de  Heer,  F.J.,  Jansen,  R.H.J.,  and  van  der  Kaay,  W.  (1979)  J.  Phys.  B ,  12:979-1002. 


where  nXe  and  are  number  densities  of  neutral  xenon  atoms  and  ambient  Ions, 

respectively,  and  V  is  the  relative  velocity  between  Xe  and  the  ambient  Ions.  Suppose  the 
ambient  Ions  have  a  zero  velocity  and  the  xenon  atoms  travel  at  a  velocity  V  equal  to  the 
maximum  vector  sum  of  the  satellite  velocity  (7.4  km/s)  plus  the  gas  exit  velocity  (0.4  km/s). 
We  obtain  V  »  7.8  km/s.  The  neutral  gas  density  is  assumed  to  be  109/cc.  The  result  of  Eq.  (2) 
is 


i\  »  3.7x102/cc  (Day) 


and 


nj  «=  l.^KP/cc  (Night)  . 


Therefore  the  ratios  p  of  seed  Ions  to  ambient  ions  are  approximately: 


p(Day)  m  3.7% 
p(Night)  =  1.2%  . 


The  estimated  results,  p(Day)  and  p(Night),  would  be  higher  If  one  takes  Into  account 
ambipolar  diffusion  and  other  mechanisms  of  electron  or  Ion  confinement,  because  these 
mechanisms  may  slow  down  the  escape  of  electrons  or  ions.  More  accurate  measurement  of 
cross  sections  would  also  help  to  Improve  the  estimated  result. 


3.4  Release  Rate 

In  thermodynamics,  the  velocity  V  of  molecules  of  mass  M  at  temperature  T  °K  is  given  by13 


vm.£dwvv).(Mi) 


1/2 


(3) 


13  Reidl,  P.C.  (1988)  Thermal  Physics,  2nd  edition,  Oxford  University  Press,  Oxford. 
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We  have  assumed  a  Maxwellian  distribution  fT(V)  to  obtain  the  result  In  Eq.  (3).  For  a  0°C 
xenon  gas,  kT  =  1.38xl0~16  x  273  gm  cm2/s2.  The  xenon  mass  M  Is  131.3x1 .67xl0"24  gm.  Using 
these  numbers  as  Input,  Eq.  (3)  gives  the  xenon  gas  velocity  V(Xe)  =  210  m/s  at  the  exit  point. 
Similarly,  for  the  C02  gas,  Its  molecular  weight  being  44,  the  gas  velocity  V(C02)  =  362  m/s. 

In  thermodynamics,  the  pressure  p(M,T)  Is  given  by 


where  p  (=  nm)  Is  the  mass  density,  and  n  Is  the  number  density  of  the  gas.  For  a  10 
atmosphere  xenon  gas,  p  =  107  gm  cm-1s-2.  The  mass  density  p  (which  Is  given  by  p  =  3P/V2)  of 
the  xenon  gas  is  p(Xe)  =  0.068  gm/cc.  Similarly,  for  the  C02  gas,  p(C02)  =  0.023  gm/cc. 

The  mass  release  rate  R,  which  Is  given  by  R  =  AVp,  depends  on  the  cross  section  A  of  the 
nozzle.  For  a  nozzle  of  1/4  cm  radius,  A  =  x(l/4)2  cm2.  Using  the  above  values  of  V,  p,  and  A  as 
Input,  we  obtain  R(Xe)  =  280  gm/s  or  84  gm/(300ms);  R(C02)  =  162  gm/s  or  48.6  gm/(300ms). 

The  particle  number  release  rates  NR  are  calculated  to  be  NR(Xe)  =  1.27xl024/s  and  NR(C02)  = 
3.8xl024/s. 

Some  Important  design  parameters  In  the  ARGOS  CIV  experiment  are  listed  in  Table  2.  The 
values  In  square  brackets  correspond  to  77°  gas  temperature. 


Table  2.  Important  Design  Parameters  for  the  CIV  Experiment  on  ARGOS 


Satellite  Velocity 

7.4  km/s 

Satellite  Altitude 

Approx  835  km 

Satellite  Orbit  Inclination 

98.7° 

Gas  Pressure  at  Nozzle 

Approx  107  gm  cm_1s'2 

Nozzle  Radius 

Approx  0.25  cm  (In  design) 

Gas  Temperature 

0°C  (77°) 

Mass  Release  Rate  (Xe) 

280  gm/s  (249  gm/s) 

Mass  Release  Rate  (C02) 

162  gm/s  [144  gm/sj 

Number  Release  Rate  (Xe) 

l. 27x10^/3  [1.14x1024/s) 

Number  Release  Rate  (C02) 

2.23x1024/s  [1.97x1024/s1 

Gas  Velocity  (Xe) 

210  m/s  [237  m/s) 

Gas  Velocity  (C02) 

362  m/s  [410  m/s) 

3.5  Release  Directions 


Since  the  spacecraft  velocity  Vs  (7.4  km/s)  is  much  higher  than  the  gas  velocity  Vg  (0.2  to 
0.4  km/s).  the  vector  sum  V  (=Va  +  Vg)  is  dominated  by  Va  (Figure  5  upper).  In  the  CIV 
experiment,  the  kinetic  energy  available  for  the  excitation  of  modified  two-stream  plasma 
instability  is  given  by  the  velocity  component  Vx  perpendicular  to  the  magnetic  field  B 
(Figure  5  lower).  For  the  CIV  process  to  occur,  Vx  has  to  exceed  the  critical  ionization  velocity 
V.  of  the  gas  released.  The  direction  of  the  earth's  magnetic  field  lines  at  Maui  is  given  by  1 1° 
declination  and  38.7°  inclination.  For  given  magnitudes  of  Va  and  Vg,  the  component  Vx 
depends  on  the  directions  of  Va  and  Vg  with  respect  to  B;  Vx  is  maximum  when  Va  and  Vg  are 
parallel  to  each  other  and  exactly  perpendicular  to  B. 


VS 

SPACECRAFT  VELOCITY 


Figure  5.  (Upper)  Vector  Diagram  of  the  Satellite  and  Gas  Velocities. 
(Lower)  Geometry  of  the  Gas  Velocity  Pitch  Angle  <t> 
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For  example,  suppose  the  satellite  is  traveling  horizontally.  Then,  Vs  subtends  an  angle  I 
with  the  magnetic  field  B,  I  being  the  magnetic  field  Inclination.  For  a  given  gas  release  angle 
0  (between  V3  and  Vg),  the  velocity  component  Vjl  and  the  velocity  pitch  angle  <j>  (between  V  and 
B)  can  be  calculated  (Figure  6).  In  this  case,  V±  exceeds  V.(Xe),  Implying  that  the  kinetic  energy 
requirement  for  CIV  Is  satisfied,  at  all  values  of  the  release  angle  8. 


3.6  Townsend's  Criterion 

Townsend  s  criterion14  states  that,  to  sustain  a  discharge,  an  electron  has  to  Ionize  a 
neutral  at  least  once  during  Its  lifetime  x  in  the  interaction  region.  In  the  literature  of  CIV, 
four  Townsend  criteria,  I  to  IV,  have  been  suggested.15  In  criterion  I,  x  Is  governed16  by  the 
escape  of  the  hot  electrons  along  the  magnetic  field  lines.  In  criterion  II,  t  is  the  contact  time 
xc  of  a  neutral  cloud17  traversing  a  magnetic  field  line.  The  longer  the  contact  time  is,  the 
more  important  Is  I  compared  to  II.  In  our  proposed  CIV  experiment,  the  duration  xc  of  the 
neutral  gas  released  is  of  the  order  of  seconds  which  is  much  longer  than  the  electron  heating 
time  xH  (order  of  ms)  obtained  by  using  18  the  estimate  xH  ■>  30/a>LH.  The  other  Townsend 
criteria  are  less  restrictive.  We  consider  criterion  I  as  follows: 

Using  a  conical  beam  model  of  half  angle  <p/2,  Lai  and  Murad16  have  formulated  the 
distance  within  which  Townsend's  criterion  I  Is  satisfied.  For  the  range  of  release  rates  In  the 
proposed  experiment,  the  estimated  distance  within  which  CIV  might  occur  is  of  the  order  of  a 
few  meters  from  the  nozzle.  Figure  7  shows  a  typical  plot  of  critical  distance  for  a  jet  of  3 
km/s,  cone  half  angle  30°.  and  molecular  weight  32  amu. 

For  the  xenon  gas  considered,  the  velocity  Is  about  7.6  km/s  In  the  space  frame,  the  cone 
half  angle  Is  30°.  and  the  release  rates  NR  have  been  estimated  In  Section  4.  Tor  simplicity,  we 
let  the  cross  section  a  of  the  electron  Impact  Ionization  of  xenon  be  o(E)  =  10“ 16  cm2  for  typical 
hot  electron  energies  E  =  10  to  30  eV  as  expected  In  a  CIV  situation.  The  Townsend  critical 
distance  X,,  =  66  m  from  the  nozzle  point.  Beyond  Xc,  the  neutral  gas  density  Is  too  low  to 
satisfy  Townsend's  criterion  I  for  CIV. 


14  Townsend,  J.S.  (1915)  Electricity  In  Gases,  Clarendon  Press,  Oxford. 

15  Lai,  S.T.  and  Murad,  E.  (1992)  Inequality  conditions  for  critical  velocity  Ionization  space 
experiments,  IEEE  Trans.  Plasma  Set.  in  press. 

16  Lai,  S.T.  and  Murad,  E.  (1989)  Critical  Ionization  velocity  experiments  In  space,  Planet. 
Space  Set.  37:7/865-872. 

17  Brenning,  N.  (1982)  Comment  on  the  Townsend  condition,  presented  at  the  Workshop  on 
Alfven's  Critical  Velocity  Effect,  Max-Planck  Instltut  fur  Extraterrestrische  Physik,  Garching, 
W.  Germany. 

18  Tanaka  M.  and  Papadopoulos,  K.  (1983)  Creation  of  high  energy  tails  by  means  of  the 
modified  two-stream  instability,  Phys.  Fluids,  26:1697. 
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GAS  RELEASE  ANGLE  8  (deg) 


Figure  6.  Perpendicular  Component  Vj_  of  Gas  Velocity  In  Space  and  the  Pitch  Angle  <)> 


We  now  derive  a  universal  formula  of  Townsend  distance  X,  for  all  masses,  assuming  that 
the  cross  sections  of  Ionization  for  each  species  are  about  the  same.  From  Eqs.  (3)  and  (4).  the 
mass  density  p  Is  given  by 


Figure  7.  Numerical  Results  for  for  Various  Hypothetical,  Conical, 
Gas  Release  Rates  In  Grams  per  300  ms  [from  Lai  and  Murad,  1989] 
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The  number  release  rate  N0  is  given  by 


N 

°  M  M 


(6) 


In  the  conical  beam  model,  the  number  density  n(x)  at  distance  x  from  the  exit  point  is  given 
by:12 


N 

n(x)  = - a- - . 

jtx2Vtan  ((p/2) 


(7) 


Townsend's  criterion  requires  that  the  electron  impact  ionization  rate  v(x)  and  the  electron 
transit  time  t(x)  should  satisfy  the  inequality: 


(8) 


where 


v(x)  =  n(x)<j(ve)ve. 


(9) 


and 


t(x)  =  2x  tan(<p/2)  u'1 , 


where  <p  is  the  full  cone  angle  and  ve  the  electron  velocity.  Using  Eqs.  (5-9),  the  universal 
Townsend  critical  distance  xc  satisfying  Eq.  (8)  is  given  by 


,  „  a  _  l-5q/lp 

c  nV  tan((p/2)  2tan(cp/2)/cT  (11) 


14 


The  formula  (Eq.  (9))  is  universal  in  the  sense  that  it  is  independent  of  the  molecular  mass  M. 
For  example,  the  critical  distance  xc  for  a  <p  =  60°  beam  is 


x  = 

C 


1.3  —  o 
kT 


(12) 


where  x,.  is  in  cm,  A  in  cm2,  p  in  gm  cm-1s-2  and  kT  in  gm.cm2/s2.  For  example,  we  take  A  = 
rr(0.25)2  cm2,  p  =  10  atm  =  107  gm  cm-1s-2,  o  =  10-16  cm2,  kT  =  273  x  1.38  x  10~16  gm-cn^/s2. 
The  Townsend  critical  distance  xc  is  66.7  m  from  the  exit  point.  This  value  of  xc  is  the  same 
for  Xe  or  C02,  or  any  gas,  as  long  as  the  ionization  cross  section  is  about  the  same. 


3.7  Altitude  Effects 


Heretofore,  all  CIV  space  experiments  were  conducted  at  around  400  to  500  km  altitude.19 
There  are  advantages  for  conducting  CIV  experiments  at  around  800  km  instead.  The  energy 
loss  to  line  excitation  of  the  ambient  neutrals20  would  be  less  at  800  km.  A  recent  model 
calculation  by  Lai  et  al21  show  that  indeed  the  energy  loss  to  line  excitation  at  800  km  is 
significantly  less  than  at  500  km  (Figure  8).  Furthermore,  the  energy  loss  to  ionization  of 
ambient  neutrals  is  also  reduced  at  800  km  where  He  (e<)>=  24.6  eV)  is  the  dominant  species. 
Since  the  ambient  neutrals,  being  thermal  and  not  beam-like,  do  not  have  a  significant  net 
velocity,  their  ionization  would  not  lead  to  any  significant  feedback  of  directed  kinetic  energy 
to  the  beam-plasma  instability  and  therefore  the  cyclic  process  of  CIV.  Both  line  excitation 
and  ambient  ionization  are  energy  loss  mechanisms. 

To  design  an  effective  CIV  space  experiment,  it  is  advantageous  to  avoid  the  above  energy 
loss  mechanisms  as  much  as  possible.  We  suggest,  therefore,  that  conducting  CIV  experiments 
at  800  km  rather  than  at  the  conventional  500  km  altitude  offers  a  better  likelihood  for  the 
CIV  experiment  to  succeed,  because  the  effect  of  both  of  the  above  energy  loss  mechanisms  is 
less  severe  at  800  km.  There  is,  however,  a  disadvantage  for  CIV  experiments  at  800  km 
compared  with  500  km;  charge  exchange,  which  is  important  for  seed  ionization,  is  less  at  800 
km  (Figure  9). 


19  Tobert.  R.B.  (1990)  Review  of  critical  velocity  experiments  in  the  ionosphere,  Adv.  Space 
Res.,  10:7/47-58. 

20  Newell,  P.T.  and  Torbert,  R.B.  (1985)  Competing  atomic  processes  in  Ba  and  Sr  injection 
critical  velocity  experiments,  Geophys.  Res.  Lett.,  12:835. 

21  Lai,  S.T..  Murad,  E.,  and  McNeil,  W.J.  (1991)  Altitude  effects  in  CIV  experiments  in  space. 
Planet.  Space  Set,  39:12  1707-1713. 
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Figure  8.  An  Example  of  Energy  Loss  to  Line  Excitation  (for  Ba  Gas)  as  a 
Function  of  Altitude  (from  Lai  et  al  1992] 
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ALTITUDE  (km) 


Figure  9.  An  Example  of  Charge  Exchange  Rate  r(h)  (for  NO  Gas)  at  Altitude  h 
Normalized  to  r(hQ)  at  h0  =  300  km  (from  Lai  et  at  1992) 


3.8  Spectral  Lines 

Electrons  energized  to  energies  above  100  eV  have  been  observed  in  previous  CIV  space 
experiments.19  If  this  happens,  then  we  expect  Xe+.  Xe+  has  a  number  of  excited  states  that 
would  be  visible  from  the  ground  by  observing  the  inter-state  emissions.  The  principal 
emitters  would  be  the  transitions  between  excited  states  of  Xe+.  For  example,  a  scenario  for 
conducting  the  experiment  is  to  release  Xe  in  sunlight,  where  the  Lyman-a  line  Is  not  enough 
to  cause  ionization.  However,  once  Xe  is  ionized  by  electron  Impact  (that  is,  CIV),  the  Lyman-a 
line  would  pump  the  Xe+  to  an  upper  state  from  which  it  would  relax  to  an  intermediate  state. 
For  example,  electron  impact  of  Xe  (by  CIV)  might  be  expected  to  give  rise  to  Xe+(2P°1/2).  This 
state  is  metastable,  with  a  lifetime  of  47  ps.  It  may  radiate  when  it  falls  to  the  ground  state 
Xe+(2p°1  l^2)  at  987  nm.  There  are  other  possibilities.  For  example,  if  the  CIV  discharge  forms 
Xe+6p(2p°i  i/2).  then  we  can  observe  the  transition  Xe+6p(2P°1  1/2  )  -»  Xe+5P5(2P°1  x/2)  at  30  nm. 
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3.9  Mass  Loading 


Haerendel22  23  pointed  out  that  the  newborn  ions  created  by  CIV  in  a  neutral  beam  can 
drag  forward  the  ambient  plasma  and  the  ambient  magnetic  field  lines.  The  forward  drag  is  a 
result  of  the  momentum  coupling  between  the  beam  ions  and  the  ambient  ions. 

To  estimate  the  critical  beam  density  nXe+  for  momentum  coupling  to  be  significant,  we 
equate  the  ion  mass  in  a  beam  pulse,  modeled  as  a  sphere  of  radius  R,  traveling  with  a  velocity 
V  and  that  in  an  ambient  flux  tube  of  length  D  (Figure  10).  The  pulse  transit  time  T  is  T  =  R/V. 
Following  Haerendel22  we  assume  that  the  length  D  is  given  by  the  distance  traveled  by  the 
Alfven  wave  in  the  same  duration  T.  That  is,  we  assume  D  =  TVA  where  VA  is  the  Alfven 
velocity.  Thus.  D  =  R(VA/V).  A  simple  estimate  of  the  momentum  balance  can  be  obtained  by 
considering  the  tube  to  be  moving  at  30  percent  of  the  velocity  V  of  the  beam: 


MXE+nXe+R  =  ManaR(VA/V)  x  30  percent 


(13) 


that  is. 


nXe+  =  na(Ma/MXe+)(VA/V)  x  0.3. 


(14) 


where  M  denotes  ion  mass,  n  ion  density,  and  a,  ambient  species.  It  is  interesting  to  note  that 
Eq.  (14)  is  independent  of  the  beam  length  R. 

Taking  typical  parameters  at  800  km  altitude,24  na  *  104/cc,  VA  «  2000  km/s,  MH/MXe+  « 
144,  V  »  7.8  km/s.  we  obtain  nXe+  »  lxl07/cc.  Thus,  for  a  neutral  Xe  gas  density  nXe  of  109/cc 
with  V  =  7.8  km/s.  a  fractional  ionization  of  a  =  nXe+/nXe  *=  1  percent  would  invoke  mass 
loading  to  such  an  extent  that  the  relative  velocity  Vr  is  reduced  to  about  Vx70  percent  =  4.5 
km/s.  Further  increase  in  ionization  a  (beyond  a  =  0.01)  may  reduce  the  relative  velocity  V, 
further,  rendering  the  crv  efficiency  of  about  1  percent  nearly  zero. 

In  the  above  example,  mass  loading  alone  already  puts  a  ceiling  on  the  maximum 
ionization  of  the  xenon  beam.  Other  energy  loss  mechanisms  lower  the  ceiling  even  further. 


22  Haerendel.  G.  (1982)  Alfven's  critical  velocity  effect  tested  in  space.  Z.  Naturforsch.. 
37A:728. 

23  Haerendel,  G.  (1983)  The  role  of  mementum  transfer  to  the  ambient  plasma  in  critical 
velocity  experiments.  Active  Expriments  in  Space,  Eur.  Space  Agency  Spec.  PubL,  ESA  SP- 
195:245. 

24  Kashe.  M.A.  (1969)  The  Ionosphere  and  Its  Interaction  with  Satellites.  Gordon  and  Breach, 

try. 
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Figure  10.  Momentum  Coupling  Between  Beam  Ions  and  the  Ambient  Ions  in  a  Flux  Tube 


3.10  Computer  Simulation  of  Xenon  CTV 


By  using  the  particle-in-cell  plasma  simulation  method25  one  can  simulate  the  explicit 
time  evolution  of  plasma  dynamics.  Earlier  works  on  computer  simulations  of  CIV  have  been 
reported  by  Machida  and  Goertz.26  Following  the  approach  of  Machida  and  Goertz,25  two  other 


25  Birdsall,  C.K.  and  Langdon,  B.G.  (1985)  Plasma  Physics  via  Computer  Simulation,  New 
York: McGraw-Hill  Book  Co. 

28  Machida,  S.  and  Goertz,  C.K.  (1986)  A  simulation  study  of  the  critical  ionization  velocity 
process.  J.  Geophys.  Res.,  91:11965. 
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CIV  simulation  codes9  27  have  been  written  independently.  These  two  codes,  using  better  cross 
sections  and  taking  into  account  full  arrays  of  colllsional  mechanisms,  yield  the  time 
development  of  CIV  as  a  result  of  the  interplay  between  collisional  and  plasma  effects. 

A  limitation  of  the  above  works  is  the  use  of  one-dimensional  space,  although  their 
velocities  are  three-dimensional.  Another  limitation  is  the  practice  of  assuming  the  ion  to 
electron  mass  ratio  to  be  100.  These  limitations  are  necessary  to  reduce  the  tremendous 
computation  time  (even  on  modem  large  scale  computers)  required  to  simulate  tens  of 
thousands  of  interacting  particles.  Despite  these  limitations,  the  simulation  results  have  been 
useful  for  providing  useful  physical  insights  into  the  detailed  dynamics9  26-27  of  CIV. 

The  purpose  of  tills  work  is  to  provide  physical  insights.  We  have  used  the  explicit  method 
as  described  in  Machida  and  Goertz26  and  McNeil  et  al9.  We  assume  that  seed  ions  are  initially 
present  in  the  xenon  cloud  (or  beam).  They  are  formed  as  a  result  of  charge  exchange  between 
Xe  neutrals  and  ambient  0+.  Before  we  present  the  CIV  plasma  simulation  results,  we  first 
present  a  Monte  Carlo  simulation  of  time  evolution  of  the  build  up  of  xenon  seed  ions  in  the 
neutral  xenon  cloud  as  a  result  of  charge  exchange  (Figure  11). 

In  our  CIV  plasma  simulation,  the  electrons  and  ions  are  represented  by  "clouds  of  charge" 
moving  through  one  another  in  response  to  the  electric  field  created  by  charge  separation, 
which  is  due  to  the  "tie  down"  of  electrons  by  the  magnetic  field  lines.  The  electrons  are  pulled 
to  some  extent  along  with  the  neutrals  through  the  polarization  drift,  but  due  to  their  short 
gyroradius  are  accelerated  mainly  along  magnetic  field  lines.  In  a  time  scale  of  less  than  a 
quarter  of  an  ion  gyroperiod,  the  ions  traverse  "straight  across"  the  field  lines  as  if 
unmagnetized.  In  this  fast  time  scale,  the  ions  accomplish  the  transfer  of  much  of  their 
kinetic  energy  to  the  plasma  instability  waves  which  in  turn  "heat"  the  electrons  rapidly. 

The  simulation  is  one-dimensional  in  displacement,  allowing  an  electric  field  in  only  one 
direction.  The  magnetic  field  is  applied  tn  a  direction  almost  perpendicular  to  the  electric 
field,  with  a  small  angle  0-  V(m/M).  where  m  and  M  are  the  electron  and  ion  masses 
respectively.  This  angle  is  chosen  for  maximum  growth  rate  of  the  modified  two-stream 
instability.  We  have  taken  into  account  the  excitation  and  ionization  of  metastable  states, 
charge  exchange,  line  excitation,  and,  of  course,  electron  impact  ionization.  The  simulation 
parameters  are  shown  in  Table  3. 


27  Person,  J.C..  Resendes,  D.,  Fetschek,  H.,  and  Hastings.  D.E.  (1990)  Effects  of  collisional 
processes  on  the  critical  velocity  hypothesis,  J.  Geophys.  Res..  95:  4039. 
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Table  3.  Simulation  Parameters 


Initial  Number  of  Particles 

1024 

Number  of  Cells 

128 

Ion  Mass/Electron  Mass 

100 

Initial  Mass  Frequency/Electron  Gyrofrequency 

1 

Tilt  Angle  of  Magnetic  Field 

6° 

Beam  Velocity/Critical  Velocity 

1.5 

Reference  Collision  Rate/coe 

0.2 

kV/coe  (Primary  Mode) 

0.04 

Ambient  Ion  Velocity 

0.0 

Typical  results  are  shown  in  Figure  12.  Electron  (upper  panel)  and  ion  (lower  panel) 
energies,  normalized  by  the  beam  energy,  are  plotted  as  functions  of  time  normalized  by  the 
initial  plasma  frequency.  The  neutral  beam  travels  in  time  from  left  to  right.  The  electrons 
are  gaining  energy  from  the  plasma  waves.  As  the  electron  energy  reaches  the  Ionization 
energy  e<i>,  Ionization  of  the  neutrals  begin  to  occur.  The  new  electrons  are  of  low  energy  when 
created  and  they  are  energized  to  higher  energies  later. 


3.11  Physical  Interpretation 

We  observe  several  interesting  points  in  the  simulation  results  shown  in  Figure  12.  (1)  The 
energy  and  population  of  the  electrons  increase  with  time.  This  is  a  manifested  symptom 
attributed  to  CIV.  (2)  The  ion  energy  spreads.  (3)  There  Is  a  new  observation  heretofore  not 
reported  in  CIV  simulations  studies.  We  observe  that  the  growth  of  Ionization  takes  off  when 
the  combined  seed  Ionization  and  electron  Impact  generated  Ionization  densities  equal  the 
ambient  plasma  density  approximately.  (4)  Taking  a  typical  nighttime  plasma  frequency  of 
1 06  Hz  at  800  km,  the  simulation  shows  that  it  would  take  about  2  ms  for  CIV  to  take  off.  For 
a  200  m/s  beam  velocity,  this  means  about  4  m  from  the  nozzle. 

The  behavior  (3)  may  be  explained  by  means  of  the  standard  linear  theory  of  beam-plasma 
instability.28  In  the  low  beam  density  regime  (o^  <  <  coa),  where  (Oj,  and  coa  are  beam  and 
ambient  plasma  frequencies  respectively,  the  modified  two-stream  instability  dispersion 
relation  yields  the  solution: 


co, 

to  =  kV-  ~1 
r  i3 


(15) 


28  Mikhailovskii,  A.B.  (1974)  Theory  of  Plasma  Instabilities,  Vol,  1,  Consultant  Bureau,  N.Y. 
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where 


(16) 


where  cdlh  Is  the  lower  hybrid  frequency:  r  and  1  denote  the  real  and  Imaginary  parts 
respectively. 


8  =  sin&^M  , 


(17) 


and 


t 


(18) 


In  the  higher  beam  density  regime  (a^  >  coa),  the  solution  is  given  by 


co 

r 


(19) 


co 


( 


Ui  ' 


(20) 


Physically,  the  ionization  rate  does  not  take  off  until  the  total  beam  plasma  density  (seed  + 
electron  impact  ionizations)  approximately  equals  the  ambient  plasma  density. 

As  a  corollary  to  the  above  theorem,  if  one  can  artificially  seed  a  neutral  beam  by  artificial 
means  (such  as  by  using  samarium,29  for  example)  to  speed  up  the  increase  of  beam  ion 
density,  it  would  trigger  CIV  more  quickly.  We  believe  that  this  point  is  worth  pursuing  for 
further  research. 


29  Lai,  S.T.,  Murad,  E.,  and  McNeil,  W.J.  (1992)  Amplification  of  critical  velocity  ionization 
by  associative  ionization,  J.  Geophys.  Res.,  97:  (No.  A4)  4099-4108. 
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4.  CARBON  DIOXIDE  GAS  RELEASE 


The  purpose  of  the  C02  gas  release  in  the  proposed  experiment  on  ARGOS  is  not  CIV  but  to 
observe  signatures  of  collisional  excitation  of  COa  with  atmospheric  O  and  H  atoms.  The 
reason  that  C02  will  not  undergo  CIV  in  this  experiment  is  that  the  critical  ionization  velocity 
V.(=  7.9  km/s)  of  C02  exceeds  the  C02  gas  velocity  in  the  space  frame.  The  magnitude  of  the 
maximum  vector  sum  of  the  satellite  velocity  Vs  (7.4  km/s)  and  the  gas  velocity  Vg  (0.4  km/s)  is 
only  I  V9+Vg  I  =  7.8  km/s.  Even  in  sunlight.  C02,  cannot  be  photoionized  by  the  Lyman-a  line  of 
sunlight  because  the  ionization  energy  [<I>(C02)  =  13.8  eV]  of  COa  exceeds  the  energy  (hv  =  10.2 
eV)  of  Lyman-a. 


4.1  Reaction  with  Atmospheric  Species 

The  800  km  altitude  of  ARGOS  is  ideal  for  studying  two  important  reactions:  the  reaction 
of  O  atoms  and  the  reaction  of  H  atoms.  In  the  first  case,  vibrationally-excited  CO  is  produced, 
while  in  the  latter  case  vibrationally-excited  OH  is  produced: 


C02  +  O  ->  CO(v)  +  02  (21) 

C02  +  H  ->  CO  +  OH(v)  (22) 


These  reactions  give  rise  to  emissions  at  the  fundamentals  of  CO(4.6p)  and  OH(2.7p)  as  well  as 
the  overtones. 


4.2  Reaction  Energetics 


The  kinetic  energy  available  for  a  collisional  reaction  (such  as  Eq.  (21)  or  Eq.  (22))  between 
two  masses  is  the  kinetic  energy  ECM  in  the  center  of  mass  frame.  By  conservation  of 
momentum. 


^space  (1+Mb/MJEcm  . 


(23) 
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Figure  13  shows  the  range,  (Eb,Ef)  of  kinetic  energy  ECM  available  for  the  reaction  [Eq.  (21)1, 
the  upper  limit,  Ef,  corresponding  to  the  forward  release  and  the  lower  one,  Eb,  to  the  backward 
release.  As  a  consequence  of  excitation  of  the  spectral  lines  in  the  energy  range,  (Eb,Ef),  the 
optical  signatures  of  the  forward  and  backward  releases  are  expected  to  be  different. 

Similarly,  Figure  14  shows  the  energy  range  for  reaction  [Eq.  (22)).  Figure  15  shows  some  of 
the  spectral  lines30  of  the  reaction  (Eq.  (21)]  in  the  range,  (Eb,  Ef). 


4.3  Monte  Carlo  Simulation 

Setayesh31  has  calculated  the  emissions  of  reactions  [Eq.  (21)  and  Eq.  (22))  and  used  a 
numerical  plume  simulation  code32  to  obtain  the  spatial  distribution  of  the  emission  in  the 
plume  of  such  a  release  (Figures  16, 17).  Setayesh's  results  indicate  that  both  of  these  processes 
[Eq.  (21),  Eq.  (22)]  would  be  observable  from  the  ground. 

Electron  impact  excitation  of  C02  may  also  occur.  In  the  C02  gas  release,  even  though  the 
relative  velocity  between  the  gas  and  the  ambient  plasma  is  just  below  the  critical  ionization 
velocity,  V..  plasma  instability  may  occur.  The  gas  can  form  an  ion  beam  in  two  ways:  (1) 
charge  exchange  between  C02  and  the  ambient  ions,  and  (2)  ambient  ion  reflection.  The 
ambient  ions  are  0+  when  in  sunlight  or  H+  at  night. 

Such  an  ion  beam,  which  is  generated  by  non-CIV  processes,  can  generate  some  CIV-like 
signatures  even  when  the  beam  velocity  is  below  the  CIV  requirement.  It  can  also  amplify  a 
CIV  process  if  the  latter  is  taking  place.29  The  ion  beam,  traversing  the  ambient  magnetized 
plasma,  can  excite  plasma  instability,  thereby  energizing  the  electrons.  Some  of  the  energetic 
electrons  may  be  able  to  excite  or  even  ionize  the  gas,  even  though  a  sustained  cyclic  ionization 
process  is  not  expected  because  of  the  subcritical  velocity  V(C02)  used  in  the  experiment. 


30  Orient,  O.J.,  Chujian,  A.,  and  Murad,  E.  (1990)  Reaction  and  electronic  excitation  in 
crossed-beams  collisions  of  low  energy  0(3P)  atoms  with  H20  and  C02,  Phys.  Rev.  Lett., 

63:No  19,  2359-2361. 

31  Setayesh.  A.  (1991)  A  parametric  study  of  the  release  of  C02  in  space.  Report  PL-TR-91- 
2052.  Phillips  Laboratory.  Hanscom  AFB,  MA  ADA-236271. 

32  Elgin,  J.B..  Cooke,  D.C.,  Tautz,  M.F.,  and  Murad.  E.  (1990)  Modeling  of  atmospherically 
induced  gas  phase  optical  contamination  from  orbiting  spacecraft,  J.  Geophys.  Res.,  95:12197- 
12208. 
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5.  CONCLUSION 


In  our  proposed  experiment,  xenon  and  carbon  dioxide  gases  will  be  released  from  the 
ARGOS  satellite.  The  xenon  gas  release  will  be  used  to  test  Alfven's  hypothesis  of  CIV.  If  CIV 
occurs,  we  expect  the  energetic  electrons  to  excite  the  Xe+,  with  subsequent  radiation.  The 
carbon  dioxide  gas  release  will  be  used  to  test  reactions  with  atmospheric  O  and  H  atoms.  The 
radiation  will  be  observed  at  Maui  and  possibly  on  the  satellite. 


COLLISION  OF  C02  +  0 


2.0  4.0 


ecm  <eV) 

Figure  13.  The  Center  of  Mass  Energy  Range,  Bracketed  by  the  Forward 
and  Backward  Release  Cases,  for  the  Collision  COz  +  O 
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COLLISION  OF  C02  +  H 


ECM(eV) 


Figure  14.  The  Center  of  Mass  Energy  Range,  Bracketed  by  the 
Forward  and  Backward  Release  Cases,  for  the  Collision  C02  +  H 
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Figure  15.  Experiment  (-)  and  Calculated  ( — )  Emission  Spectra 
in  the  1B2->1A1  Transition  of  COa  for  the  Collision  C02  +  O 
[from  Orient  et  al,  1990] 
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Figure  16.  Results  of  Monte  Carlo  Simulation  of  C02  +  O  CO*  +  02  Showing 
the  CO  Vibrational  Emission  at  4.6  nm 
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re  17.  Results  of  Monte  Carlo  Simulation  of  C02  +  H  ->  OH*  +  CO  Showing 
)H  Vibrational  Emission  at  2.7  |im 
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